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TITLE OF T HE INVENTTOM 

BACKGROTmm ni? nrrs E mWNTJr™ 

The invention relates to seafloor electromagnetic (EM) surveying for oil and 
other hydrocarbon reserves. 

Geophysical methods for mapping subterranean resistivity variations by 
various forms of EM surveying have been in use for many years [,, 2, 3, .0]. In these 
methods, eieotrie field detectors are piaced on the seafloor a, carefittiy chosen 
posttions a. rar,ges up ,o about ,0 km from an electiomagnetic source. Detector 
s.gr.is measured at the detectors are sensitive to variations in subterraneart strata 
oonfignratton resistivity beneath the area being surveyed. However, EM surveying 
was not widely though, of as a technique tha, couid be app.ied to finding hydrocarbon 
reservoirs. 

More recently, it was proposed to use EM surveying to find hydrocarbon 
res=rvo,rs..An early proposal by Statoil was to use me vertical current flow 
components to detect hydrocarbon layers [4, 5], since i, is these components that am 
sensmve to the presence of a thin resistive layer. Inis was based on the understanding 
that a subterranean stiata configuration that inchtdea a resistive hydrocarbon layer 
embedded within .ess msistive sediments wil. give rise to a measurable enhancement 
of the e.entiic field amplitude compared to a subterranean stra* confiscation 
compnstng only water-bearing sediments. Tne Statoil pmposal was to collect data 
fiom detector iocations which are in-line with (i.e. end-on to, me axis of a horizontal 
eleemc dtpole (BED) antenna so that the ga,vanic mode, tha, should be most sensitive 
,0 the presence of a buried high resistivity ,ayer, domina,es. However, i, was 



- ~» two . xauwcver, n was j 

satisfied ma, the Statoi! method could not provide retiabie result, since the in-line 



j 0 , , . . wumaw j oxnv,c uic m-nne 

dati .collected , mcapable of distinguishirtg b«ween a thin buried hydrocarbon layer 
of htgh resistivity siaated in less res . st . ve ^ ^ ^ one hand, and a non- 
hydrocarbon bearing rock formation in which me strata exhibits increasing resistivity 



I 



with depth, on the other hand, the latter being a common feature of many large scale 
sedimentary structures. 

It was then proposed to use the EM surveying method according to Sinha [12] 
for finding hydrocarbon reservoirs [9, 13] and it was then confirmed that this method 
works well in practice for finding hydrocarbon reservoirs [6, 7]. The essence of the 
Sinha method is to normalise the in-line data with equivalent data for the same 
source-detector pair locations collected in an orthogonal geometry where the inductive 
mode dominates the response, referred to as a broadside geometry. In the broadside 
geometry the axis of the HED dipole antenna of the source is perpendicular to a line 
between the detector and source. EM surveying of a hydrocarbon reservoir applying 
the Sinha method is now described in more detail. 

Survey data is collected by using a surface vessel to tow a submersible vehicle 
carrying a HED antenna over a survey area. The HED antenna broadcasts a source 
electromagnetic signal into the seawater. Detectors are located on the seafloor over the 
survey area and measure a signal in response to EM fields induced by the HED 
antenna. The amplitude of the detector signals is sensitive to resistivity variations in 
the underlying strata configuration and this is used to determine the nature of the 
subsea structure. In order to successfully map subterranean resistivity variations, the 
orientation of the current flows induced by the source electromagnetic signal must be 
considered [6]. The response of seawater and subterranean strata to the source 
electromagnetic signal is different for horizontally and vertically flowing induced 
current components. For horizontally flowing current components, the coupling 
between the layers comprising the subterranean strata is largely inductive. This means 
the presence of a thin resistive layers (which is indicative of a hydrocarbon layer) does 
not significantly affect the detector signal at the seafloor since the large scale current 
flow pattern is not affected by the thin resistive layer. On the other hand, for vertical 
current flow components, the coupling between layers is largely galvanic (i.e. due to 
the direct transfer of charge). In these cases even a thin resistive layer strongly affects 
the detector signals at the seafloor since the large scale current flow pattern is 
interrupted by the resistive layer. 
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" comp °" ents of induced cumnt flow whi <" - — 

mm*, to the presence of a dun resiaive layer> sofc ^ 
or eating a hyd^carbon ^ „ „ ^ ^ J - 

s,gna,s a, the detectors arising ^ to presence . ^ ^ 

&om ^ effecK wMch ^ from o(her reaiist;c - 

^de e ae response of , he subKnanean state to horizonta[ 

coup,ed) and vesica, (i.e. vertical* coup,ed) induced current flows M 

10 and . , An . de0,rOmagne,iC "«» "* - « HED antenna generates both inductive 
and galvanrc current flow modes, with the rela tive strengti, of each mo de depending 
—e-detector geometry. At detector iocations which are broadside «o I HED 
«-» drpo e axis, the inductive mode dominates the response. At detector iocaZ 

S, 9, ,0]. Accordrngiy, the response of the subterranean strata to vertica, indue d 

" ^ ^ ^ — * — - — orientation I I 
detector and the response of the subterranean strata ,o horizonta, induced cur™, 

detector. Data from both geometric configurations is required 

bn ■ H theref ° re imP<>tta,,, fai8ning 3 Practica) EM -rv=y &r detecting 

sitr: layers mtos ^ ^ to — - ^ ~ 

detector confutations „ which the coupiing between iayers is ,ar g e,y inductive due 
o ^ _ (in which case the survey has tittie sensitivity to tiT^ ^ 
dun resrstive iayer, and Arose in which tire coupihtg between layers isTJI 
.alvantc due to vertica, entente (in which case blocking of tire passage „ fth ^ 

zi::c ,Mds ,o a - y wwch is ~* — * - ~~ ' 

Figure , shows in plan view an examp.e survey geotnetiy accruing to tire 
Stnha metitod. TTrere are sixteen detectors 2S, and tirese are !aid I in a 
on a section of seafloor 6 above a subterranean hyd^carbon reservou T " 
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hydrocarbon reservoir 56 has a boundary indicated by a heavy line 58. The orientation 
of the hydrocarbon reservoir is indicated by the cardinal compass points (marked N, 
E, S and W for North, East, South and West respectively) indicated in the upper right 
of the figure. To perform a survey, a source such as a HED antenna, starts from 
location 'A' and is towed along a path indicated by the broken line 60 through 
location 'B' until it reaches location «C\ which marks the end of the survey path. As 
is evident, the tow path first covers four parallel paths aligned with the North-South 
direction to drive over the four "columns" of the detectors. This portion of the survey 
path moves from location 'A' to location 'B'. Starting from location l B\ the survey 
path then covers four paths aligned with the East- West direction which drive over the 
four "rows" of detectors. Each detector is thus passed over in two orthogonal 
directions. The survey is completed when the source reaches the location marked £ C\ 
During the towing process, each of the detectors 25 presents several different 
. orientation geometries with respect to the source. For example, when the source is 
directly above the detector position Dl and on the North-South aligned section of the 
tow path, the detectors at positions D5, D6 and D7 are at different ranges in an end-on 
position, the detectors at positions D2, D3 and D4 are at different ranges in a 
broadside position and the detector at positions D8 and D9 are midway between. 
However, when the source later passes over the detector position Dl when on the 
East-West aligned section of the tow path, the detectors at positions D5, D6 and D7 
are now in a broadside position, and the detectors at position D2, D3 and D4 are in an 
end-on position. Thus, in the course of a survey, and in conjunction with the 
positional information of the source, data from the detectors can be used to provide 
details of the source electromagnetic signal transmission through the subterranean 
strata for a range of distances and orientations between source and detector. Each 
orientation provides varying galvanic and inductive contributions to the signal 
propagation. In this way the continuous towing of the source can provide a survey 
which samples over the extent of the subterranean reservoir. 

• The Sinha method has been demonstrated to provide good results in practice. 
However, it has some limitations. 



# 



Firstly, since the two modes cannot be easily separated there will generally be 
a level of cross-talk between then, a, a detector and this can lead to ambiguities in the 
results 
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Secondly, in order to obtain survey data from bod, in-line and broadside 
5 geometrres, the HED antenna needs to present two different orientations a. each 
source location. Tins requires the surface vessel „ make 
broadcast locations and can lead to long and complex towpath patterns 

Thirdly, me survey can only provide the best data possible a, discrete source 
locations. This „ because rf fte geometr . c requiremMts ^ a ^ ^ 

whteh motate mat, a, any point during me survey, data can only be optimally collected 
from those detectors to which the HED antenna is arranged either in-line or broadside 
At other orientations, separation of the inductively and galvanically coupled signals 
becomes more difficult, and resulting data are less reliable. For instance, referring to 
tire figure, when me HED antenna is a, a point on the tow path directiy above the 
detector marked Dl and on the North-South aHgned section of the tow path, in-line 
data can only be coUected from me detectors marked D5, D6 and D7, whilst broadside 
data can on,y be collected torn the detectors marked D2, D3 and D4. The other 
detectors (for examp.e mose marked D8, D9 and D 10) provide only marginally ussM 
mfcrmanon at mis point of tire survey because of the comp.ex mixing of the 
galvamcaUy and inductively coup,ed modes. Furthermore, if, for example, the HED 
antenna is a, tire location identified by reference numeral 57 in the figure, which is on 
a North-South aligned section of the tow path, in-line data cat, be collected from the 
detectors marked D3, D8, D9 and D.0, bu, broadside data cannot be coUected from 
any of tire detectors. Since both broadside and Wine dafa are required f or optimaI 
analysrs, the best data possib.e with the sc.ua. detector array shown in the figure can 
only be collected from pointa along the tow path where the source is directly above 
one of the detector locations. 

In summary, with the Sinha method, the time during which good quality data 
can be eofieoted represents omy a small fiaction of the overalt time taken to perform a 
survey. W™, h addition to tire aurvey being time-inefficient, it is necessary to 
accurately follow a complex tow path which has to complement the detector layl 



and tire detector, Selves mus, a!so be c^My accurateIy ^ 
d.fficul„es in controHtag both the position and «be orientation of a towed source 
antenna, coupied with this need to accurately foUow a particuiar tow path relative to 
tire detector grid, is one of tbe major sources of error in surveys of these kind. The 
disadvantages associated with are survey constraints imposed by the Sinha method are 
the pnce to pay for resolving the ambiguities inherent in the Statoil method 
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meti, d f ^ inVemi ° n *~ iS Pr0VMed -e y 

m«hod for surveying „ M , hat . s ftought m . ^ (o contefa a 

hydros on reservoir, comprising: ^ . souroe electromagnetic 

location; dating . ^ sigMl . a ^ locaUon fa * 
and ob,armng survey data indicative of phase difference between firs, and second 

>0 By comparing phase measurements of different components of tine de.ee.or 

?* MPara,i0n — * CM be d «~^ which is sensitive .0 ft. presence 

of a hydrocarbon iayer or reservoir within a subterranean stia* configuration. The 
or no, a phase separation anontaly, and hence me presence or no. of a 
h^ocar on layer, can be determined with a singie scarce orientation. Ttae is no 

r« ,s wi * b,own methods tased °- &r ** » * 

nrr : ource orien,ations - - * — 

qutckiy and wtthon, needing ,0 accurateiy con.ro, the soume orientation. Furthetmore 
^cauae of titis msensMv* of a phase measnremen, ,0 me reiative source orientatioT 
citable da. nonaction is no, limned ,0 specific source iocation and 
20 geometrtes, as ta me case when coHecting broadside amplitude da,a, and a 

ZTT1 is also possible ,o use 311 «* <"•*». <* .0 

ha ve first, second and third components. 

bv obse^ ^ SeC ° nd "* ^ ^ direCti ° nS ^ Preferab1 ^ ^onai, since 
by obs rvtng geometiicaliy independen, component of me de,ec,or signais, thete is 

presence of a hydrocarbon reservoir is accordingly increased. 
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The source e,ecttomagnetic signa, can be broadcast from an antenna mounted 
on a submersible vehicle, or from a static iocation, such as within a borehole, or from 
an oil ,or gas platform. 

5 obta- ^ d ~ 8ne,iC «« b = « different fluencies to 

obtant survey dato „ a pluraIity of d . fferent fcquendes Morco ^ ^ 

O m TTZT 1 °" te emitted ^ 3 ™^ rf P^rab.y between 

0.0, Hz and ,0 Hz. The method can be advantageously repeated over the same survey 
area ustng different frequencies of source electromagnetic signal. Lower frequencies 
are generally preferred. By probing the subterranean strata a, a number of differen 
frequences of soume electromagnetic signal, i, is possible to obtain improved vertical 
resolution of structures within the subterranean stiata configuration 

The source signs! can be from a horizontal electric dipole. Such a signa, can be 
« using eating equipment, and also al,ows rel a«ve,y simp, inversion 

; The invention also provides a method of analysing result from an 

e ^magnetic survey of an area that is thought or known to contain a subterranean 
hydrocarbon reservoir, comprising: providing survey data indicative of phase 
difference between firs, and second components of the detector signal solved lug 
fc, and second duections respective,* exacting the phase differences from the 
survey dam; and determining a metric from the phase differences ,ha, is predictive of 
the presence or absence of hydrocarbon. " * 

The phase differences can be extracted by rotational* transforming the survey 
data from an mstrument frame to a source frame. 

The invention a!so provides a computer pmgram product bearing machine 
re able tnstiucuons for dementing the anaiysis method, and a computer loaded 
with such a computer program product. 

The invention finther provides a method of fanning an electromagnetic 
survey of an ama tha, is though, or known to contain a subterranean hydrocarbon 
~, comprising: creating a mode, of the area to be surveyed including a seafloor, 
rock fonnation containing a posttuated hydrocarbon reservoir beneath the aeafloor 
and a body of water above the seafloor; setting va,ues for depth be.ow the seafloor of 
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the postulated hydrocarbon reservoir and resistivity structure of the rock formation- 
and performing a simulation of an electromagnetic survey in the model to obtain from 
the model phase differences between first and second components of a detector signal 
resolved along first and second directions respectively. 
5 Repeated simulations for a number of source-to-detector distances and 

frequencies can be performed in order to allow optimum surveying conditions in 
terms of source-to-detector distance and frequency of EM signal for probing the 
hydrocarbon reservoir to be selected when performing an electromagnetic survey The 
effects of differing detectors array configurations and source tow paths can also be 
10 modelled. 

The invention also provides a computer program product bearing machine 
readable instructions for implementing the planning method, and a computer loaded 
with such a computer program product. 

15 
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BRIEF DESCRTPTTQN OF THF nPAWTwrc 

For a better understanding of the invention and to show how the same may be 
carried into effect reference is now made by way of example to the accompanying 
5 drawings in which: 

Figure 1 is a schematic plan view showing an example survey geometry 
following prior art principles in which sixteen detectors are laid out on a section of 
seafloor above a subterranean reservoir; 

Figure 2A shows in schematic vertical section a surface vessel undertaking an 
10 EM survey; 

Figure 2B is a plan view detailing a polar coordinate system; 
Figure 3 shows in schematic vertical section a model uniform background 
subterranean strata configuration; 

Figure 4 shows in schematic vertical section a model hydrocarbon-layer 
1 5 subterranean strata configuration; 

Figure 5A shows a graph plotting calculations of phases of different 
components of detector signals seen during a model electromagnetic survey of the 
subterranean strata configurations shown in Figures 3 and 4; 

HgUre 5B * ^Ph P'omng differences in the phases shovm in Figure 

5A; 

Figure 6 shows in schematic vertical section a model of a non-hydrocarbon 
bearing subterranean strata configuration; 

Figure 7 shows a graph plotting calculations of phases of different components 
of detector signals seen during a model electromagnetic survey of the subterranean 
strata configurations shown in Figures 3, 4 and 7; 

Figure 8 shows a graph plotting differences in the phases shown in Figure T 
Figure 9A is a schematic plan view showing an example survey geometry 
according to an embodiment of the present invention in which sixteen detectors are 
laid out on a section of seafloor above a subterranean reservoir; 



25 
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Figure 9B compares the signal coverage of the prior art method and the 
method of the invention; 
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Figure 10 shows a graph plotting calculations of phases of different 
com pon e nts of detector signals seen during a model electromagnetic survey of the 
subterranean strata configurations shown in Figures 4 at two different source 
electromagnetic signal frequencies; 

5 Figures 1 1 A and I IB show graphs p.otting differs in the phases shmvn in 

rigure 10; 

Figure 12A shows a graph plotting emulations of differences in phases 
between different components of detector signals seen during a mode, electromagnetic 
survey „f several hydrocarbon-layer subterranean stiata configurations with an 
1 electromagnetic source signal frequency of 0.25 Hz; 

Figure 12B shows a graph plotting calculations of differences in phases 
between different components of detector signals seen during a mode, electromagnetic 
survey of several hydrocarhon-layer subterranean strata configurations with an 
electromagnetic source signal frequency of 0.5 Hz; 

Figure 12C shows a graph plotting calculations of differences in phases 
between different component* of detector signals seen during a model electromagnetic 
survey of several hydrocarbon-iayer subterranean strata configurations with an 
electromagnetic source signal frequency of 1.0 Hz; 

F ' 8 "" 120 shows * Plotting calculations of differences in phases 

between dtfferent components of de.ec.or signals seen during a mode, electromagnetic 
survey of severe. hydrocarbon-.ayer subterranean sttata configurations win, an 
electromagnetic source signal frequency of 2.0 Hz; 

Figure 13A shows a graph plotting calculations of phases of different 
components of detector signals seen during a model electromagnetic survey of several 
uniform background subterranean strata configurations; 

Figure 13B shows a graph plotting differences in the phases shown in Figure 

Figure ,4 shows a graph plotting calculations of maximum observed 
dtfferences in phases between differen. component of de,ec,or signal seen during a 
model electromagnetic survey of a hydrocarbon-layer subterranean strata 
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configuration, as a function of hydrocarbon-layer resistivity, and for several 
electromagnetic source frequencies; 

Figure 15A shows a graph plotting calculations of phases of different 
components of detector signals seen during a model electromagnetic survey of the 
subterranean strata configurations shown in Figures 3 and 4; 

Figure 15B shows a graph plotting differences in the phases shown in Figure 

i. jA j 

Figure 16A shows a graph plotting calculations of phases of different 
components of detector signals seen during a mode! electromagnetic survey of the 
subtetranean strata configurations shown in Figures 3 and 4; and 

Figure 16B shows a graph plotting differences in the phases shown in Figure 

1 0/\ ? 
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A method of electromagnetic surveying for oil and other hydrocarbon reserves 
. descnbed which doe, no. require separate data acquisition of the response of a 
subterranean strata configuration to inductive* and galvanicafly coupied modes He 
new method can be performed using pre-existing survey equipment 

Figure 2A schematicaily shows a surface vessel .4 undertaking EM surveying 
of a subtenanean stiate configuration in a way that is suitable for coUecting survey 
data for carrying o„« the invention. Th. subterranean stiate configuration includes an 
overburden ,ayer 8, an underburden ,ayer 9 and a hydrocarbon ,ayer (or reservoir) a 
The surface vesse, ,4 floats on me surface 2 of the seawater 4. A deep-towed 
submerge vehic.e ,9 carrying a HED antenna 2, is attached to the surface vessei ,4 
by an umbilical cable ,6 providing an electrical, optical and mechanical connection 
between the de=p-,owed submersible vehicle 19 and the surface vessel 14 The HED 
antenna broadcasts a source electromagnetic signal into the seawater 4 

One or more remote detectors 25 are located on the seafloor 6. Each detector 
25 tndu es an instrument packages 26, a detector antenna 24, a floatation device 28 
and a bafias, weigh, (no, shown). The detector antenna 24 measures a detector signal 
nr response ,o EM fields induced by the HED antenna in the vicinfty of me detecter 

25 fte amphtrde of the detector signais is sensitive ,o resistivity variations in me 
underiymg strata conflguratjon ^ ^ ^ ^ ^ 

for ater analysis. The detector antenna 24 in .his exampie comprises ,wo 
Orthogonal d.pole antennae arranged .o detec firs, and second component of ft. 
electttc field in a horizontal plane, i.e. one which is paralle. ,„ the seafloor 6 

The de.ec.ors record two (or three) orthogonal component of «he seafloor 
electee fie,d as raw date. The raw survey d*a are men analysed, after recovery of me 
de,ec,ors and transfer of the raw date in,o a suftabte compter. Initially a spectia! 
analysts ,s performed to remove tite componen, of me signs, which corresponds ,o 
source transmission, as is conventional. The survey date are then combined wfth 
source and receiver navigation data, again as is conventional. Then me survey date are 
processed ,o rotete the electric fieids from an ■inatinmenf ftame (i.e. component 
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paralle. ,o to receiver dipoles of to detector) ,o to •source' feme (i.e. radia! and 
aztmutol components referenced «o to source-receiver geometry). This is a new 
processing step specific to the present invention. 

Figure 2B is a schematic plan view detailing a polar-coordinate system which 
ts nsed ,o describe to princip.es of to new method. The origin of to coordinate 
system is positioned at to centre of to HED antenna shown in Figure 2A, and zero- 
azrmuth is aligned parailel to to dipole axis of to HED antenna, as indicated in 
Ftgnre 2B (to HED antenna in this Figure no, drawn to scale). In Figure 2B, a single 
detector 25 is shown positioned at a range of R km from to origin, and a, an azimuth 
of * . The orthogonal dipo.e antennae comprising to detector antenna 24 are 
arbitrary oriented in to horizontal plane as indicated in to figure 

Since to phase difference from to source of electromagnetic radiation from a 
honzontal dipo.e source is an azimuthally symmcteic subterranean strata configuration 
ts largely independent of azimuth O, phase measurements recorded at to detector 
antenna 24 are .argely independent of to azimtrtoi position of to detector 25 shown 
« Figure 2B. This al.ows phase data to be coltected e q ual,y over a wider ran«e of 
source-detector orientations man is possib.e with amp.itude data, and any inaccuracies 
m to measurement of to azimuthal position of to detector in to coordinate system 
shown in Figure 2B have a lesser effect 

' In to following examples, to two components of detected electric field also 

known as detector signal, for which to phase is measured axe a radial component and 
an aztmutol component. The radial component is to, component of to electric field 
resolved atong a direction parallel to a line connecting to source location and to 
detector iocation, and marked E, in Figure 2B. The azimuto. component is to, 
component of to electric field resolved a,ong a direction perpendicular ,„ a line 
connecting to source tocation and to detector location and in a horizontal plane, and 
marked E. in Figure 2B. The components of to detected electric field atong these 
directions is determined from to angular orientation of to ortogonal dipole 
antennae comprising to detector antenna 24 relative to to line joining to source 
location and to detector location. This can be easi.y determined using standard 
tnstrumemation, such as, for example, active or.passive sonar to determine to relative 
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positions of .he source .ocation and «. detector locatioDj ^ . 
determine the detector antenna orientation. 

to order to show how the respective phases of two spatia, components (e g 
radta. p and azimuthal * components) of me electric field can be used to detect the 
presence of a subterranean hydrocarbon reservoir, nunterica, forward modeffing of tire 
tad descttbed by Chave and Cox t „J is applied to ^ mode , ^.^^ 
strata configurations. 

Figure 3 shows in schematic vertical section a mode, background subtetranean 
sttata configuration. The configuration comprises a section of seafloor .06 beneath a 
km depth of seawater 104. The seawater has a resistivity of 0.31 fin, Beneath the 
seafloor 106 is a unifonn half-space sedimentary structure with a resistivity of 1 flm 
tire low resistivity being primarily due to aqueous saturation of pore spacas This' 
background subterranean stiata configuration extends uniformly downwards for an 
• mftotte extent. Also indicated in Figu, 3 are a HED antenna 2,, and a detector 25 
such as mose shown in Figure 2A. The distance between tire HED antenna and the' 
detector (i.e. the range) is R km. The azimuthal position of me detector Native to the 
onentation of the HED antenna is arbitral due to the insensitivity of the phase 
component of the detected electric field signals to azimuth. 

Figure 4 shows in schematic vertical section a model hydrocarbon-.ayer 
subterranean strata configuration. A section of seafloor .06 lies beneath a 10 km 
depth of seawater ,04 which has a resistivity of 0.3, flm. The strata configuration 
beneath the seafloor ,06 comprises a , km thick overburden layer ,08, representing 
sedtments, arranged above a hydrocarbon layer 1 ,2. The overburden !ayer 108 has a 
reststivity of , nm, again, primarily due to aqueous saturation of pore spaces. The 
hydrocarbon .ayer 1 ,2 is 0.1 km thick, and has a resistivity of ,00 fim. The re,a«ve,y 
htgh resistivity of the hydtocarbon layer is due to the presence of non-conducting 
hydrocarbon within pore spaces. Below the hydmcarbon layer 1,2 is a sedimentary 
underburden ,ayer ,09, which, as for the overburden layer, has a resistivity of , Om 
The underburden ,ayer extends downwardly for an effectively infinite extent 
Accordingly, except for the presence or absence of the hydrocarbon layer „2 the 
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background subterranean St rata configuration of Figure 3 and the hydrocarbon-layer 
sub.errar.ean strata configuration of Figure 4 are identica,. A HED antenna 2, and a 
detector 25 are again shown as in Figure 3 . 

5 azimuth^T ^ ^ * P ' 0ninS m ° deUed «*- 9 » f - 

—ai component of dae de,ec te d elecrnc field for born tire backgmund 

ub.erranean stote configuration and dae hydrocarbon-iayer subrerranean stra.a 

-mgurauon nrode.s shown in Figures 3 and 4 respectively as a Sanction of range R 

^ Phase , measnred relative .o a source ^magnetic signal .ransmided b y dae 

of 0.5 Hz. Tb. rad,a, and azimuma, co mp onen.s of dae de.ec.ed elecric fleid for dae 
baokground subrerranean strata configuration are marked e p ° and 8.= respectively and 
dae corresponding component of the de te c.ed elecric field for dae hydrocarbon-layer 

show daa. 6p s, ... epR ^ a „ ^ ^ fa ^ * 
.5 Howeve, „ „ ^ clear «ha t dae ra,e of phase advanoe is less for dae hydrocarbon-' 
layer subterranean sti-a* oonf.gm.don daan for the background subterranean strata 

ToT* T ^ ° f baCk8IDUnd SUb '~ — dae 
Phase of both dae radial and azhnutiaa, component advances a. a ra«e of around W 

per to. In addition, a. ranges beyond abou. 2 km, dae azimuma, component ^ 
0 on^enfly lags the radial component 0/ by mmd 25 , _ ^ 

subter^nean stra* configuration, however, me behaviour is somewhat differ!, 
Beyond around 5 km, me azunmha, componen. C again lags the radial componen. 

1 1 0 ^ ^ ^ ° f ^ C °~ - * - °^und 

on y 10 per km. This is significantiy lower .han that seen wim dae background 

subterranean stea configmation. Furthermore, a. ranges between around 2 km and 5 

km, the difference in phaae hereon azimuma, componen. 9.° and me radia, 

componen. e p varies significant*. A phase separation anomaly is seen which varies 

d t; : '° °' ^ differe ° Ce ^ "» '^ — " -ponen B of 
detected electric field to a maximum of almost 60\ 

Kgnre 5B shows a graph plotting, f or both me radial and azimuthal 
components, the difference in chase f) R n" k«,„. t , 

in pnase u -e between the hydrocarbon-layer and 
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background subterranean strata configurations as a function of range R. The difference 
m the radial components is marked e p R - 0p B , and the difference in azimuthal 
components is marked 9^. The different rates of phase advancement seen with 
the hydrocarbon-layer and background subterranean strata configurations is apparent 
m the negative gradient of the curves beyond around 3 km. The differing behaviour at 
and ranges (between around 2 km and 5 km) is apparent from the separation of the 
curves over this range. 

These differences in phase behaviour, namely the relatively slow advancement 
a Phase of both radial and azimuthal components when a reservoir is present, and the 
strong range-limited variation in phase between the radial and azimuthal components 
seen at mid-ranges when the reservoir is present, provide two useful characteristics 
with which to determine the presence or absence of a hydrocarbon layer within an 
otherwise uniform background. 

For the practical application of controlled source electromagnetic surveying to 
hydrocarbon exploration, it is necessary that other common subterranean strata 
configurations do not lead to a behaviour similar to that seen in the hydrocarbon-layer 
subterranean strata configuration model. In particular, it is important to be able to 
. distinguish between subterranean strata configurations which include a thin 
hydrocarbon layer and non-hydrocarbon containing subterranean strata configurations 
that have increasing resistivity with depth. 

Figure 6 shows in vertical section a highly schematic model of a non- 
hydrocarbon containing subterranean strata configuration. This subterranean strata 
configuration exhibits increasing resistivity with depth, which is a common feature of 
many large scale sedimentary structures. Due, for example, to increasing expulsion of 
conducting seawater with depth from rising overburden pressure. As with the 
background and hydrocarbon-layer subterranean strata configurations described 
above, in the non-hydrocarbon bearing subterranean strata configuration a section of 
seafloor 106 lies beneath a 10 km depth of seawater 104. The strata beneath the 
seafloor 106 comprise a series of sedimentary layers of increasing resistivity. A first 
layer 1 10 has a unifonn resistivity of 1 Q m and a thickness of 1 km. A second layer 
1 13 has a uniform resistivity of 5 Qm and a thickness of 1 km. A third layer 1 14 has a 
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uniform resistivity of 50 Qm and a thickness of 1 km. Beneath the third layer 1 14 is a 
fourth layer 116 which has a resistivity of 100 Qm and extends downwardly for an 
infinite extent. A HED antenna 21 and a detector 25 are again shown as in Figure 3. 

Figure 7 shows a graph which is similar to and will be understood from the 
5 description of Figure 5 A above, but which also includes modelled curves determined 
for the non-hydrocarbon bearing subterranean strata configuration. The modelled 
curves marked 9 P B , 9«> B , 0 p * and 8** are the same as those shown in Figure 5A. The 
curves marked 9 P S and 9<, s show the radial and azimuthal components of the detected 
electromagnetic field seen with the non-hydrocarbon bearing subterranean strata 
10 configuration. 

The behaviour of the variation in phase of the detected electric field as a 
function of range beyond around 5 km is broadly similar for the hydrocarbon-layer 
subterranean strata configuration and the non-hydrocarbon bearing subterranean strata 
configuration. There are moderate differences in gradient between the models for the 
1 5 example shown, and in some circumstances this may allow the two subterranean strata 
configurations to be distinguished. (In other examples, the gradients are almost the 
same.) However, even if measurable, the value of the gradient is likely to be a fairly 
unreliable indicator of subterranean strata configuration in practice. This is because 
different absolute values of resistivity, for instance a more or less resistive 
10 hydrocarbon layer in the hydrocarbon-layer subterranean strata configuration, or a 
more rapidly increasing resistivity with depth in a non-hydrocarbon bearing 
subterranean strata configuration, are likely to lead to changes in the observed 
gradients and cause confusion between the two models. 

However, at ranges between around 2 km and 5 km, there is nothing in the 
5 phase of the detected electric fields in response to the non-hydrocarbon bearing 
subterranean strata configuration which resembles the phase separation anomaly seen 
with the hydrocarbon-layer subterranean strata configuration. The phase behaviour 
seen with the non-hydrocarbon bearing subterranean strata configuration much more 
closely resembles that of the background subterranean strata configuration across this 
) range. Accordingly, it is the phase separation anomaly (and not the gradient) which 
provides the most appropriate indicator of subterranean strata configuration. 
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Figure 8 shows a graph ploning «ha difference to toe phase AO between the 
rad,al and azimutha! components of the detected electric field for the three model 
subterranean shara configurations descrihed above as Sanction of ranged Recurve 
marked A6 in Fig^e , represents rhe difference between me curves markeu ^ and 
P m F.gure 7 (with negative values corresponding to the azhnutoal component 
laggmg me radial component). The curves marked A6* and A9* in Figure 8 
co^pondinglv represent rhe differences between ,he curves marked * and ^ and 
andGp in Figure 7 respectively. 

10 sua, ^ Ph3Se SePMa,i0n ^ ^ 11,6 yer subterranean 

suata configurahon (curve marked A8» to Figure 8) is apparent as a clear hou»h 

centred a, a range of around 3.5 km. The magnitude of me phase separation anomaly 
a. to pom, ,s almost 60". This difference to phase between the radia, and azimutoa. 
components of toe delected electric field is about 30- more negative man the largest 
drfference seen with either toe background or- toe non-hydrocarbon bearing 
ubterranean sua. configuration, Wito current technology, a phaae differences 
between the radial and azimutoa, components of toe delected electric fie,d of 10° can 
be clear* resolved. According toe presence or not of a .rough simUar to mat seen in 
F.gure 8 » easily detectoble, and able to distinguish between a hydrocarbon-layer 
subterranean stoala. configuration of toe type shown in Figure 3, and toe mode, 
0 subterranean strata configurations shown in Figures 4 and 6. 

By distributing a linear array of delectors a.ong a section of seafloor, and a. 
each one recording suitable raw data in response to a sourca electromagnetic signa, 
broadest by a horizon*, electromagnetic dipoie source, pl o« such as toose shown in 
_ F.gure 8 can be generated from toe phase informal obtained from toe raw data The 
> results of toese p, ots can then be used to indicate toe type subterranean strata 
configurahon beneato a hue joining toe source and toe detectors. Unlike previous 
survey methods, this can be done using a single dipole source and witoou, toe need 
to collect multiple data sets corresponding to different orientations of toe source 

By drstributing a planar array of detectors on a section of seafloor, and a, each 
one recording raw data in response ,o a source e.ectromagnefic signa, broadcast by a 
honzomal etoctromagnetic dipole source, plots snch as those shown in Figure 8 can be 
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generated for a number of different directions onee the phase information has been 
extracted from the raw data. Becanse of the insensitivhy of phase ,o me azimuthal 
posmon of a detector with respect to the source dipole axis, the pica a.ong each of the 
Afferent directions achievable with a planar array of detectors can be obtained 
sunulteneously, irrespective of the dipole source orientation. This allows a thorough 
two- or three-dimenaional survey ,„ be performed without even having to move the 
source. This contrasts to previous methods where a relatively long and complicated 
tow of the dipole source is required to utU.se all of the detectors in a planar array and 
then only with relatively low spatia! sampUng. White all of me detectors in a pianar 
' array can be utihsed without moving the source, in practical surveys emptying the 
new method, it is „ ke ly that the source will nonemdess be moved, such as shown in 
Ftgure 2A. Each new source position provides an entire set of useful source-detector 
geometr.es, artd so provides more comprehensive sampling of the subterranean strata 
configuration for a given number of detector, In addition, by moving me source 
surveys can be fully performed where the detectors are deployed over an area with a 
characteristic scale farger than me range of distances over which phase measurements 
can be reliably used to indicate the presence of a hydrocarbon layer. 

Figure 9A shows in plan view an area of seafloor 6 to be surveyed and which 
ts stmilar ,o drat shown in the prior arf Figure 1. There are sixteen detectors 25 for 
recordmg the phase components described above. The defector are laid out in a square 
gnd above a subterranean reservoir 56. Other detector distributions could be used 
msfead, such as other grid shapes, or distributions mat are no, in a simple grid (The 
constramts on detector placemen, patterns imposed by me amplitirde-based memods 
of bom Stefoil and Sinha are merefore lifted.) The subterranean reservoir 56 has a 
boundary indicated by a heavy ,ine 58. The orientation of me subterranean reservoir is 
mdrcated by the cardinal compass points (marked N, E, S and W for North East 
South and Wes, respectively) indicated in me upper righ, of me figure. To perform a 
survey using an embodimen, of me new method, a source star* from location 'A' and 
■a towed along a pam indicated by me broken line 120 ,o location 'B>, which marks 
the end of the survey pam. A, most points along tire tow pam, useflU da,a M be 
collected from all of me detecters. For example, when me source antenna is a, me 
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location marked by the reference numeral 57 in Figure 9A, all sixteen of the detectors 
25 are able to collect reliable data. This contrasts to the correspondingly similar 
locate shown in Figure 1, again marked by the reference numeral 57, at which point 
no useful data can be collected using previous methods. Accordingly, the tow path 
5 shown m Figure 9A, which is startlingly simple compared to that shown in Figure 1 
actually provides a much greater amount of valid data. As noted above, with the ne J 
method, rt a only necessary to know the relative positions of the source and detectors 
and the orientation of each detector antenna such that the radial and azimuthal 
components of the detected electric field can be geometrically resolved. Since the 
10 onentation of the antenna is not critical, there is no need for the tow path 120 shown 
in Figure 9A to closely follow a pattern defined by the grid of detectors. In fact it is 
preferable for the tow path to not align too closely with the north-south and east-west 
based detector grid, since for detectors in an end-on position (i.e. at an azimuth of 0' 
. m the coordinate system shown in Figure 2B), the amplitude of the azimuthal 
5 component of the detected electric field will be small for a dipole source, and the 
Phase of this component more difficult to accurately establish. In source-detector 
onentations where either the radial or azimuthal components of the detected electric 
field are small, other components of the detected electric field may be employed, for 
instance as described further below. 

3 Figure 9B is a graphical representation comparing the signal spread of the new 

Phase method to the old inline/broadside amplitude method. The example reservoir 56 
bounded by the perimeter 58 is shown. The dipole source is at an arbitrary location 57 
wxthm the reservoir with the HED antenna axis aligned W-E. In the old method good 
quahty inline amplitude data is only collectable within a narrow angular range 64 
UKhcated by W-E dark shading in the figure, and good quality broadside amplitude 
data is only collectable within a narrow angular range 62 indicated by the N-S dark 
shadxng in the figure. The angular ranges 62 and 64 need to be narrow to ensure that 
one of the inductive and galvanic signal components dominates over the other Data 
collected by detectors in the four main quadrants 66 is essentially bad data to be 
rejected from analysis. On the other hand, in the new method, the situation is reversed 
The broad quadrants 66 become the regions over which good quality data is collected 



-22- 



10 



15 



20 



25 



30 



since they are the regions in which phase can be fully decomposed into the radial and 
azimuthal signal components needed for the phase-difference anomaly measurement 
whereas the dark areas 62, 64 are angular areas where the collected data becomes 
unreliable since the magnitude of one of the radial and azimuthal signal components is 
likely to become too small causing signal-to-noise problems. 

The modelled phase responses shown in Figures 5A, 5B, 7 and 8 were all 
calculated for a horizontal electromagnetic dipole source transmitting a source 
electromagnetic signal at a frequency of 0.5 Hz. 

Figure 10 shows a graph plotting the modelled phase 8 of the radial and 
azimuthal components of the detected electric field for the hydrocarbon-layer 
subterranean strata configuration model shown in Figures 3 as a function of R for two 
different frequencies of source electromagnetic signal. The modelled radial and 
azimuthal components of the detected electric field seen in response to a dipole source 
transmitting at a frequency of 2 Hz are marked 6 P 2 » 2 and ft, 2 " 2 respectively, and the 
modelled radial and azimuthal components of the detected electric field seen in 
response to a dipole source transmitting at a frequency of 0.25 Hz are marked 0 p °- 25H2 
and 9o°- 25fe respectively. These curves, and also comparison with the modelled radial 
and azimuthal components of the detected electric field seen in response to a dipole 
source transmitting at a frequency of 0.5 Hz, marked 6 P R and ft/ in Figure 5A 
indicate a frequency dependence to the characteristics of the phase separation anomaly 
indicative of a buried hydrocarbon layer. Towards higher frequencies, the phases of 
the rad,al and azimuthal components advance faster than at lower frequencies, and the 
scale over which the phase separation anomaly characteristic of a hydrocarbon layer's 
presence occurs is also seen to be frequency dependent. 

Figure 1 1 A shows a graph plotting the difference in the phase AG between the 
radial and azimuthal components of the detected electric field for the reservoir and 
background model subterranean strata configurations described above as function of 
range R, hv response to a dipole source transmitting at a frequency of 2 Hz. The curve ' 
marked A9 R in Figure 11 A represents the difference between the curves marked ft,** 
and 9 P 2HZ in Figure 10 (with positive values corresponding to the radial component 
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lagging the azimuthal component). The curve marked A9 B represents the 
corresponding data for the background model subterranean strata configuration. 

Figure 11B is shows a graph plotting the difference in the phase A9 between 
the radial and azimuthal components of the detected electric field for the reservoir and 
> background model subterranean strata configurations described above as function of 
range R, in response to a dipole source transmitting at a frequency of 0.25 Hz. The 
curve marked A9*in Figure 1 IB represents the difference between the curves marked 
0^ Hz and 9 p 025Hz in Figure 10 (with positive values corresponding to the radial 
component lagging the azimuthal component). The curve marked A0 B represents the 
corresponding data for the background model subterranean strata configuration. 

It is clear from Figures 8, 1 1 A and 1 IB, that the range over which the phase 
separation anomaly occurs is smaller at higher frequencies. At 2 Hz (see Figure 1 1 A), 
the phase separation anomaly is centred at a range of around 3 km and occurs over a 
characteristic range of about 1 km. At 0.5 Hz (see Figure 8), the phase separation 
anomaly is centred at a range of around 3 km and occurs over a characteristic range of 
about 2 km. At 0,25 Hz (see Figure 1 IB), the phase separation anomaly is centred at a 
range of around 3 km and occurs over a characteristic range of about 3 km. For all 
frequencies, the maximum phase separation seen with the hydrocarbon-layer 
subterranean strata configuration is about 30° more negative than the phase difference 
that would be seen if the hydrocarbon layer were not present. This indicates that the 
presence of a hydrocarbon layer can be detected using a range of frequencies, each of 
which acts a probe of the subterranean strata configuration operating over a slightly 
different spatial scale. 

Figure 12A is a graph showing the effect of differing overburden thicknesses. 
The graph plots the difference in the phase AG between the radial and azimuthal 
components of the detected electric field for several hydrocarbon-layer subterranean 
strata configurations with different overburden thicknesses as function of range R. In 
this example, the source electromagnetic signal is at a frequency 0.5 Hz. Curves are 
plotted for different hydrocarbon-layer subterranean strata configurations which, while 
otherwise similar to the hydrocarbon-layer subterranean strata configuration shown in 
Figure 4, have overburden thicknesses of 0.25 km, 0.5 km, 1.0 km, 1.5km and 2.5 km. 
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The curves corresponding to each different overburden thickness are correspondingly 
marked in the figure. The curve marked 1.0 in Figure 12A is identical to the curve 
marked A8 R in Figure 8 since the overburden thickness in the model shown in Figure 
3 (and used for the modelling shown in Figure 8) is 1.0 km. The curve marked A8 B in 
Figure 12A is similar to and will be understood from the similarly marked curve in 
Figure 8. For each of the curves corresponding to different overburden thicknesses, 
the magnitude of the phase separation anomaly is roughly similar, varying from about 
55° with an overburden thickness of 0.25 km to about 65° with an overburden 
thickness of 2.5 km. Accordingly, the method is equally able to detect a thin 
hydrocarbon layer at a range of depths beneath the sea floor. It is also apparent that the 
range at which the phase separation anomaly is maximum increases with increasing 
overburden thickness. This sensitivity of the range of maximum phase separation 
anomaly to overburden thickness can allow the depth of a reservoir to be determined 
with appropriate inversion modelling and suitable data coverage. 

Figures 12B, 12C and 12D are similar to and will be understood from Figure 
12A. However, Figures 12B, 12C and 12D show the response of different overburden 
thicknesses to different frequencies of source electromagnetic signal. Figure 12B 
shows the response to a source electromagnetic signal at a frequency 0.25 Hz, Figure 
12C shows the response to a source electromagnetic signal at a frequency 1.0 Hz, 
Figure 12B shows the response to a source electromagnetic signal at a frequency 2 Hz. 
It can be seen that the phase separation anomaly is detectable with a range of 
frequencies over a range of overburden thicknesses. The magnitude of the phase 
separation anomaly is broadly similar at each of the different frequencies shown. As 
seen previously, the range over which the phase separation is apparent narrows with 
increasing frequency. 

Figure 13 A is a graph showing the effect of differing background resistivity. 
The graph plots the modelled phase 9 of the radial and azimuthal components of the 
detected electric field for a background subterranean strata configuration similar to 
that shown in Figure 3, but with different resistivity values for the uniform 
subterranean strata, as a function of range R. In this example, the source 
electromagnetic signal is at a frequency 0.5 Hz. The phase of the radial and azimuthal 
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components of the detected electric field are calculated for resistivity values of 1 fim, 
5 Om, 15 £2m, 50 D.m and 200 Sim, as marked in the figure. For each resistivity value 
the phase of the radial component of the detected electromagnetic field is shown as a 
dashed line, and the phase of the azimuthal component is shown as a solid line. The 
pair of curves corresponding to the 1 Clm resistivity value are identical to the curves 
marked 9 P B and e<p B in Figure 5A. It is clear that the resistivity value for the 
subterranean strata in a uniform model containing no hydrocarbon reservoir has a 
significant effect on the detected phase. The rate of advancement of phase with range, 
for both the radial and azimuthal components of the detected electric field, falls with 
increasing background resistivity. For example, at a resistivity of around 15 Qm, the 
rate of phase advancement with range is about 15° per km. This is similar to the rate 
of phase advancement seen with the reservoir subterranean strata configuration model 
and plotted in Figure 5A for ranges beyond around 5 km. This again demonstrates 
how absolute values of phase for each of the radial and azimuthal components can be 
an unreliable indicator of the likely presence of a hydrocarbon layer within an 
otherwise uniform resistivity background. 

Figure 13B shows a graph plotting the difference in phase A9 between the 
radial and azimuthal components of the detected electric field for each of the different 
resistivity value background subterranean strata configurations shown in Figure 13 A. 
Each curve is appropriately marked according to the resistivity value of the model to 
which it corresponds. The curve marked 1 Qm is identical to the curve marked A9 B in 
Figure 8. While the characteristic difference in phase for the radial and azimuthal 
components is depends on the resistivity of the uniform subterranean strata, none of 
the curves shown in Figure 13B display a range limited phase separation anomaly 
which, as seen in Figure 8, is indicative of the presence of a buried hydrocarbon layer. 
It is the resistivity contrast between a buried hydrocarbon layer and an otherwise 
uniform background which gives rise to the phase separation anomaly. Accordingly, 
by forming the phase difference between the radial and azimuthal components of the 
detected electric field in the manner described above, a hydrocarbon-layer containing 
subterranean strata configuration remains clearly distinguishable from a range of 
uniform subterranean strata configurations of differing resistivities. 
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Figure 14 is a graph showing the effect of hydrocarbon-layer resistivities. The 
graph plots the largest difference in phase max(A9) between the radial and azimuthal 
components of the detected electric field as a function of differing hydrocarbon-layer 
resistivity P, in a hydrocarbon-layer subterranean strata configuration which is 
> otherwise similar to that shown in Figure 3. Curves are shown for source 
electromagnetic signal frequencies of 0.25 Hz, 0.5 Hz, 1 Hz and 2 Hz, as indicated in 
the figure. For example, with a source electromagnetic signal frequency of 0.5 Hz, and 
a hydrocarbon-layer resistivity of 100 Qm, the largest difference in phase between the 
radial and azimuthal components of the detected electric field is about -58'. This 
particular value corresponds to the minimum seen in Figure 8 for the curve marked 
A0 R . Typical hydrocarbon-layer resistivities are between a few tens of Qm and a few 
hundreds of Hm. It can be seen from Figure 14 that for all typical values of 
hydrocarbon-layer resistivity, a difference in phase between the radial and azimuthal 
components of the detected electric field of at least 30' is seen for all source 
electromagnetic signal frequencies. At lower frequencies, it is even higher. 

This demonstrates that the above described method is able to detect 
hydrocarbon layers with different resistivities, and using a range of source 
electromagnetic signal frequencies. 

It can also be seen from Figure 14 that the maximum phase difference between 
the radial and azimuthal components of the detected field is greatest for a hydrocarbon 
layer with a resistivity of about 50 Qm. At resistivities above and below this value, a 
decreasing maximum phase difference is seen. 

It has thus been demonstrated that the presence of a hydrocarbon layer in a 
subterranean strata configuration can be detected by observing the phase difference 
between radial and azimuthal components of detected electric field in response to a 
source electromagnetic signal from a horizontal electric dipole source. This has been 
shown to work over a wide range of source frequencies, for differing depths of burial 
of a hydrocarbon layer and for different subterranean strata configuration resistivity 
values. 

Alternative Embodimgnte 
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Whilst in the above examples the radial and azimuthal components have been 
considered, similar techniques can also be employed using different components of 
the detected electromagnetic field. For instance, if the detectors 25 shown in Figure 2 
5 were configured to also record the phase of the vertical component of the detected 
electric field (i.e. perpendicular to both the radial and azimuthal directions), the 
vertical component could be used in combination with another component to probe a 
subterranean strata configuration. 

Figure 15A shows a graph plotting the modelled phase 8 of the vertical and 
10 azimuthal components of the detected electric field for both the background 
subterranean strata configuration and the hydrocarbon layer subterranean strata 
configuration models shown in Figures 3 and 4 respectively as a function of R. The 
phase is measured relative to a source electromagnetic signal transmitted by the HED 
antenna. In this example, the source electromagnetic signal is at a frequency 0.5 Hz. 
The vertical and azimuthal components of the detected electric field for the 
background subterranean strata configuration are marked 0 Z B and a." respectively and 
the corresponding components of the detected electric field for the hydrocarbon-layer 
subterranean strata configuration are marked 8 2 R and 9** respectively. Except for 
showing the vertical rather than the radial components of the detected electric field, 
20 this figure directly corresponds to Figure 5A. 

Figure 15B shows a graph plotting the difference in the phase AG between the 
vertical and azimuthal components of the detected electric field for the two model 
subterranean strata, configurations included in Figure 15 A. The curve marked A9 B in 
Figure 15B represents the difference between the curves marked 0 O B and 9 2 B in Figure 
15A (with negative values corresponding to the vertical component lagging the 
azimuthal component). The curve marked A9 R Figure 15B correspondingly represents 
the difference between the curves marked Q 0 R and 9 2 R in Figure 15 A. 

Figure 16A shows a graph plotting the modelled phase 9 of the vertical and 
radial components of the detected electric field for both the background subterranean 
strata configuration and the hydrocarbon-layer subterranean strata configuration 
models shown in Figures 3 and 4 respectively. The phase is measured relative to a 
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source electromagnetic signal transmitted by the HED antenna. In this example, the 
source electromagnetic signal is at a frequency 0.5 Hz. The vertical and radial 
components of the detected electric field for the background subterranean strata 
configuration are marked Q» and 9 P B respectively and the corresponding components 
of the detected electric field for the hydrocarbon-layer subterranean strata 
configuration are marked 9 Z R and 9 P R respectively. Except for showing the vertical 
rather than the azimuthal components of the detected electric field, this figure directly 
corresponds to Figure 5A. 

Figure 16B shows a graph plotting the difference in the phase AG between the 
vertical and radial components of the detected electric field for the two model 
subterranean strata configurations included in Figure 16A. The curve marked A9 B in 
Figure 16B represents the difference between the curves marked 9 P B and 9 Z B in Figure 
16A (with negative values corresponding to the vertical component lagging the radial 
component). The curve marked A9 R Figure 18B correspondingly represents the 
difference between the curves marked 6 P R and 9 2 R in Figure 1 8A. 

Figures 15B and 16B both indicate that the difference in phase between the 
vertical component of detected electric field and either of the azimuthal or radial 
components is also sensitive to the presence of a hydrocarbon layer in an otherwise 
uniform background subterranean strata configuration. The phase difference seen 
between the vertical and azimuthal components displays both a negative and a 
positive lobe compared to the background subterranean strata configuration with a 
crossover at a range of around 4 km. This would be particularly useful indicator for 
use in survey areas where background resistivity is poorly constrained. The qualitative 
behaviour of the phase of the vertical component of the detected electric field is 
approximately similar to that of the azimuthal component. However, a larger phase 
separation is seen when comparing the radial component with the azimuthal 
component than when comparing the radial component with the vertical component. 
Accordingly, the azimuthal component will generally be preferred when comparison is 
made with the radial component, unless, for instance, the magnitude of the azimuthal 
component is small, for example, where a detector is very close to an end-on 
orientation. 
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In the above description, and in Figures 5 A, 7, 10, 13 A, 15A and 15B, the 
absolute phase of various components of the detected electric field has been 
considered relative to the source electromagnetic signal phase. However, in practice, 
' since it is the relative between different components of the detected electric field seen 
at the detector which is indicative of the presence of a hydrocarbon layer, the detected 
components may be directly compared without reference to the phase of the source 
electromagnetic signal. 

Finally it will be understood that the invention is equally applicable to 
surveying of freshwater, for example large lakes, so that references to seafloor, 
. seawater etc. should not be regarded as limiting. 



Summary 

It has been demonstrated how a phase separation anomaly occurs in response 
to a hydrocarbon layer which is not seen with a background subterranean strata 
configuration. This allows the detection of subterranean hydrocarbon reservoirs and 
hydrocarbon bearing layers. The technique has many advantages over previous 
methods, for example: 

Previous techniques based on comparison of amplitude measurements require 
the collection of both end-on and broadside data for each receiver to be reliable. This 
requires multiple orthogonal survey tow paths (see Figure 1). Using a technique such 
as described above, a survey may be completed more thoroughly with a much shorter 
and less complex tow path (see Figure 9A). 

• The orthogonal towpaths required by previous methods lead to sampling of 
different parts of a target structure. Because only single source-receiver pairs are 
required for a phase-based detection of the reservoir, there is reduced interpretational 
ambiguity arising from the dimensionality of the target structure. This overcomes the 
limitation of the prior art Sinha method in which in-line data from a given receiver 
will come from one source position and the corresponding broadside data at the same 
range will generally come from a different source location. This means that the 
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structure sampled between the source and detector for the compared in-line and 
broadside data will not be the same. With the new method, this problem does not 
arise, since all the data is collected from a single source position, so both phase 
components in the processed signal are derived from sampling the same structure. 

• The above described technique is almost independent of detector azimuth 
relative to a source's dipole axis. Since the method is less dependent on the 
orientation of the source, geometry-related errors are much reduced. In order to 
decompose the detector signals into radial and azimuthal (or whichever components 
are desired) it is only necessary to know the relative positions of the source and 
detectors, and the orientation of the detector antenna. These can be easily determined 
using existing technology. 

• The phase separation seen above is range limited and can be controlled by 
varying the frequency of the electromagnetic source. If the source were to broadcast at 
several discrete frequencies (either by employing multiple source antenna or a tuned 
source for example) improved vertical resolution can be achieved. 

• For a particular source dipole transmission frequency, the range dependence of 
phase separation can be used to indicate the depth to the resistive layer. 

• Phase data are relatively insensitive to structures which are local to the 
receiver. 
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1. An electromagnetic survey method for surveying an area that is thought or is 
known to contain a subterranean hydrocarbon reservoir, comprising: 

transmitting a source electromagnetic signal from a source location; 
detecting a detector signal at a detector location in response thereto; and 
obtaining survey data indicative of phase difference between first and second 
components of the detector signal resolved along first and second directions 
respectively. 

2. The survey method of claim 1, wherein the first and second components are 
radial and azimuthal. 



3. The survey method of claim 1, wherein the first and second components are 
vertical and azimuthal.' 



4. The survey method of claim 1, wherein the first and second components are 
vertical and radial. 



5. The survey method of claim 1, further comprising obtaining survey data 
indicative of phase of a third component of the detector signal resolved along a third 
direction orthogonal to the first and second directions. 

6. The survey method of claim 5, wherein the first and second and third 
components are vertical, radial and azimuthal. 



/• The survey method of any one of the preceding claims, wherein the first and 
second directions are orthogonal. 
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8. The survey method of any one of claims 1 to 7, wherein the source 
electromagnetic signal is broadcast from an antenna mounted on a submersible vehicle 
which is towed over the survey area to move the source location. 

5 9. The survey method of any one of claims 1 to 7, wherein the source location is 
fixed. 

10. The survey method of any one of the preceding claims, wherein the source 
electromagnetic signal emits is emitted at different frequencies to obtain survey data at 

1 0 a plurality of different frequencies. 

11. The survey method of any one of the preceding claims, wherein the wherein 
the source electromagnetic signal is emitted at a frequency of between 0.01 Hz and 10 
Hz. 

15 

12. A method of analysing results from an electromagnetic survey of an area that 
is thought or known to contain a subterranean hydrocarbon reservoir, comprising: 

providing survey data indicative of phase difference between first and second 
components of the detector signal resolved along first and second directions 
20 respectively. 

extracting the phase differences from the survey data; and 
determining a metric from the phase differences that is predictive of the 
presence or absence of hydrocarbon. 

25 13. The analysis method of claim 12, wherein the first and second components are 
radial and azimuthal. 

14, The analysis method of claim 12, wherein the first and second components are 
vertical and azimuthal. 
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15. The analysis method of claim 12, wherein the first and second components are 
vertical and radial. 

16. The analysis method of claim 12, further comprising obtaining survey data 
indicative of phase of a third component of the detector signal resolved along a third 
direction orthogonal to the first and second directions. 

17. The analysis method of claim 16, wherein the first and second and third 
components are vertical, radial and azimuthal. 

18. The analysis method of any one of claims 12 to 17, wherein the first and 
second directions are orthogonal. 

19. The analysis method of claim 18, wherein the phase differences are extracted 
by rotationally transforming the survey data from an instrument frame to a source 
frame. 

20. A computer program product bearing machine readable instructions for 
implementing the method of any one of claims 12 to 19. 

21. A method of planning an electromagnetic survey of an area that is thought or 
known to contain a subterranean hydrocarbon reservoir, comprising: 

creating a model of the area to be surveyed including a seafloor, a rock 
formation containing a postulated hydrocarbon reservoir beneath the seafloor, and a 
body of water above the seafloor; 

setting values for depth below the seafloor of the postulated hydrocarbon 
reservoir and resistivity structure of the rock formation; and 

performing a simulation of an electromagnetic survey in the model to obtain 
from the model phase differences between first and second components of a detector 
signal resolved along first and second directions respectively. 
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22. The planning method of claim 21, wherein the first and second components 
are two of radial, vertical and azimuthal. 

23. The planning method of claim 21 or 22, further comprising: 

5 repeating the simulation for a number of source-to-detector distances and 

frequencies in order to select optimum surveying conditions in terms of source-to- 
detector distance for probing the hydrocarbon reservoir. 

24. A computer program product bearing machine readable instructions for 
1 0 implementing the planning method of claim 2 1 , 22 or 23 . 

25. An electromagnetic survey method substantially as hereinbefore described 
with reference to Figures 2 to 16 of the accompanying drawings. 

15 26. A method of analysing results from an electromagnetic survey substantially as 
hereinbefore described with reference to Figures 2 to 16 of the accompanying 
drawings. 

27. A method of planning an electromagnetic survey substantially as hereinbefore 
20 described with reference to Figures 2 to 16 of the accompanying drawings. 
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ABSTRACT 

ELECTROMAGNETIC SURVEYING FOR HYDROCARBON RESERVOIRS 

An electromagnetic survey method for surveying an area that potentially contains a 
5 subterranean hydrocarbon reservoir. The method comprises detecting a detector signal 
in response to a source electromagnetic signal, resolving the detector signal along at 
least two orthogonal directions, and comparing phase measurements of the detector 
signal resolved along these directions to look for a phase separation anomaly 
indicative of the presence of a buried hydrocarbon layer. The invention also relates to 
10 planning a survey using this method, and to analysis of survey data taken using this 
survey method. The first and second data sets may be obtained concurrently with a 
single horizontal electric dipole source antenna. The method is also largely 
independent of a source-detector pair's relative orientation and so provides for good 
spatial coverage and easy-to-perform surveying. 

15 
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